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Abstract 

Quasi-stationary distributions, as discussed by Darrocli & Seneta 
(1965), have been used in biology to describe the steady state be- 
haviour of population models which, while eventually certain to become 
extinct, nevertheless maintain an apparent stochastic equilibrium for 
long periods. These distributions have some drawbacks: they need not 
exist, nor be unique, and their calculation can present problems. In this 
paper, we give biologically plausible conditions under which the quasi- 
stationary distribution is unique, and can be closely approximated by 
distributions that are simple to compute. 
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1 Introduction 

The logistic growth model of Verhulst (1838) was the first to describe math- 
ematically the evolution of a population to a non-zero equilibrium, contrast- 
ing with the Malthusian law of exponential growth. Its stochastic version, a 
Markov chain X in continuous time in which X{t) represents the number of 
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individuals at time t in a population in a prescribed area A, has transition 
rates 



transitions 



q-limits (1.2 



qse 



(1.1) 



Qij = 



bi, 



_i = di + ei'^/A, 



i > 1; 
otherwise, 



where b and d are the per capita rates of birth and natural mortality, and 
there is an additional per capita death rate ex, due to crowding, at pop- 
ulation density x = i/A. The stochastic model has the drawback that its 
equilibrium distribution assigns probability 1 to the state zero, population 
extinction, irrespective of the initial state. This apparently negates the 
most valuable property of Verhulst's model, its ability to allow an equilib- 
rium other than extinction. However, if 6 > d and A is large, the popula- 
tion density X{t)/A can be expected to remain near the 'carrying capacity' 
K := [b — d)/e for a very long time, in an apparent (and often biologically 
relevant) non-extinct stochastic equilibrium. 

Darroch and Seneta (1965), building on the work of Yaglom (1947) in the 
context of branching processes, introduced the concept of a quasi-stationary 
distribution, in an attempt to reconcile these at first sight inconsistent prop- 
erties of the model. In a discrete time Markov chain X consisting of an 
absorbing state together with a single finite transient aperiodic class C, 
the limiting conditional probabilities 



qj := hm P,[X(t) = j \ X{t) G C], i,j & C, 



exist, and are the same for each i G C. The qj, j G C, also determine a 
quasi-stationary distribution, in the sense that 



(1.3) 



qk 



j€C kec 



where P := (pjk) denotes the one step probability transition matrix. If, 
however, C is countably infinite, the situation is very much less satisfactory; 
there may be no quasi-stationary distribution, or exactly one, or infinitely 
many, and determining which of these is the case may be a difficult problem. 
Even when there is a unique quasi-stationary distribution, its calculation 
can pose substantial problems, unless the equations (jl.3p happen to have an 
obvious solution, because the probabilistic definition (|1.2|) involves condi- 
tioning on an event which, in the limit as t ^ oo, has probability zero. This 
appears to make the quasi-stationary distribution unsatisfactory for typical 
biological applications. 
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In this paper, we give conditions, simply expressed in terms of the prop- 
erties of the process X, under which things are in fact much simpler. Under 
the conditions of Theorem 12. 11 there is exactly one quasi-stationary distribu- 
tion, and it can be approximated to a specified accuracy by the equilibrium 
distribution tt^ of a 'returned process' X^. What is more, under slightly 
more stringent conditions, the distribution of X{t) is shown in Theorem 12. 51 
to be close to the quasi-stationary distribution for long periods of time. 

The returned process, introduced by Bartlett (1960, pp. 24-25) and used 
by Ewens (1963, 1964) in a population genetical setting, is a Markov process 
that evolves exactly like X, up to the time at which is reached, but is 
then instantly returned to a random state in C, chosen according to the 
probability measure n. The mapping fi i-^ n'^, studied in the paper of Ferrari 
et al. (1995), is contractive under our conditions, and iterating the mapping 
leads to the unique quasi-stationary distribution m on C, which satisfies 
m = tt"^. In many practical applications, including the stochastic logistic 
model of (jl.ip when A is large, iteration is unnecessary, inasmuch as any 
distribution tt^ is extremely close to m. Furthermore, since ir^ is a genuine 
equilibrium distribution, its computation does not involve conditioning on 
sets of vanishing probability, and is hence typically much simpler. 

The main results. Theorems 12.11 and 12. 5( are proved in Section [2j In 
Section [3l as an illustration, we discuss the application of the theorems to 
birth and death processes, of which the stochastic logistic model p.ip is an 
example. These processes have the advantage of having been widely studied, 
because of their relatively simple structure, allowing our results to be easily 
interpreted; however, the theorems are equally applicable to processes with 
more complicated structure. 



Let X be a stable, conservative and non-explosive pure jump Markov process 
on a countable state space, consisting of a single transient class C together 
with a cemetery state 0. For any probability distribution /x on C, define the 
modified process X^ with state space C to have exactly the same behaviour 
as X while in C, but, on reaching 0, to be instantly returned to C according 
to the distribution fi. Thus, if Q denotes the infinitesimal matrix associated 
with X, and that belonging to X'^, we have 



2 A general approximation 



return-measure 




for i,j G C. 
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tau-def 



U-def 



equilibrium-bnd (2.4) 



In this section, under a rather simple set of conditions, we show that the 
quasi-stationary distribution m of X is unique, and can be approximated in 
total variation to a prescribed accuracy by the stationary distribution of X'^, 
for an arbitrary choice of /U. We give a bound on the total variation distance 
between m and vr^ that is expressed solely in terms of hitting probabilities 
and mean hitting times for the process X, and which is the same for all fi. 
The bound is such that it can be expected to be small in circumstances in 
which the process X typically spends a long time in C in apparent equilib- 
rium, before being absorbed in as a result of an 'exceptional' event. If the 
bound is not, as it stands, small enough for practical use, it can be improved 
geometrically fast by iteration of the return mapping ^ tt^. 
Our basic conditions are as follows. 

Condition A. There exist s £ C , p > and T < oo such that, uniformly 
for all k e C, 

(i) pk := Pfc[X hits s before 0] > p; 

(ii) Efe[r{s,o}] < T < oo. 

Here, and refer to the distribution of X conditional on X{0) = k, and 

(2.2) TA ■■= mf{t > 0: X{t) G A, X{s) ^ A for some s < t}, 

the infimum over the empty set being taken to be oo. Condition A (i) can be 
expected to be satisfied in reasonable generality; Condition A (ii), although 
satisfied by the stochastic logistic model, is not so immediately natural. 
We now introduce the quantity 



(2.3) 



u ■■= ^qko/{qA{r{k,o})}- 



To interpret the meaning of U, observe that a renewal argument for X^, 
with renewal epochs the visits to any specific j G C, shows that 



TT^U) < l/{g,E,(r{,-o})}. 



In particular, if X has a quasi-stationary distribution m, it follows from 
that 

U > ^vr™'(i)gio = '^'m{i)qiQ = Xm, 

where is the rate at which the A-process, starting in the quasi-stationary 
distribution m, leaves C: Pm[A(t) G C] = e~^"'^. Thus U acts as a com- 
putable upper bound for any A^- Note that p,T and U are all quantities 
that can reasonably be bounded using a knowledge of the process X. 
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quasi-approx 



s-mean 



In the remainder of this section, we show that the quasi-stationary dis- 
tribution m exists, is close to any vr^, and well describes the long time 
behaviour of X prior to absorption in 0, as long as UT/p is small enough. 
Our first main result is the following. 

Theorem 2.1 Suppose that Condition A is satisfied, and that 2UT/p < 1. 
Then X has a unique quasi- stationary distribution m, and, for any proba- 
bility measure fx on C, we have 

d^y{m,'!r^') < 2UT/p. 

Remark. Of course, for the theorem to imply that vr^ is a sharp approxima- 
tion to m, one needs U to be small enough (and therefore certainly finite). 
In many applications, X can only jump to from a small number of states 
in C, and, if the quasi-equilibrium really behaves like a genuine equilibrium 
for long periods of time, the quantity IEfc(r{fc o})) for each such k, can be 
expected to contain a large contribution from paths that, after leaving k, 
spend a very long time 'in equilibrium' in other states of C, before either 
returning to k or being absorbed in 0. In such applications, as in the next 
section, these two features combine to make U small. 

To prove the theorem, we first need some preparatory results. We first 
show that, under Condition A, the mean time to hitting the state s is uni- 
formly bounded, for all return processes X'^, and for all initial states. 

Lemma 2.2 Under Condition A, for all probability measures fx on C and 
for all r & C, we have 

^rrj^.y < T/p < (X, 

where is defined similarly to ta, but with the process Xi^ in place of X. 
Proof. Recursively define 

^{5,0} — ^{s,oy 

:= inf{t > rfjfo^^: X'^(t) € {s,0}, ^ {s,0} for some rfjfo^^ <u<t], 

J>2, 

and, for j > 1, let Zj := /[X^(r[^' = 0, 1 < I < j], taking Zq = 1. Then 
it follows that 
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s-second 



Now Ej.Ti'jQi < T, by Condition A(ii), and, for j > 2, 

by Condition A(ii), where T denotes the u-field of events up to the 
stopping time T^fo^ ■ Then, for j > 1, 



E{Zj\T^,,,-i} < {l-p)Zj_i, 

{s,0} 



by Condition A(i). Hence, for j > 1, and for any r G C, it fohows that 
and so ^rT^gj < T/p, as required. I 

It follows in particular from Lemma 12.21 that KgT^^y < T/p < oo, so 
that X^^ is positive recurrent on C; denote its stationary distribution by vr^. 
Then, for any /: C — )• M bounded, set 

/•oo 

h%3) ■■= - {E,/(Xn-vr^(/)}dt, ieC. 
To show that the integral is well-defined, note that 

/•oo 

\h%j)\ < / 2\\f\\d^^{C{Xnt)\XnO)=jUndt, 



where || • || denotes the supremum norm. The latter integral is finite provided 
that Es{(r|'_,j)^} < oo, by the coupling inequality (Lindvah 2002, (2.8)) and 
from Pitman (1974, Corollary 1, (1.23) with r = 2). That this is the case 
follows from the next lemma. 

Lemma 2.3 Under Condition A, for all probability measures /i on C , we 
have 

E,{(rf^p2| < ^_ 
Proof. Writing r := r|*^p note that 



(poo \ poo / poo \ 

/ I[T>t]dt] =2 I[T>t]i I[t >u]du] dt. 
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Now, from Lemma 12.21 and by the Markov property, we have 
E (^j^ I[t > u] du ^ < {T/p)I[t > t], 

where J-j^ denotes the history of X^^ up to time t. Hence, taking expectations 
in (j2.5p . it follows that 

IEs(r2) < 2(T/p)E,T < 2{T/pf, 

again from Lemma 12.21 completing the proof. I 



It is shown in the proof of Theorem 12. 51 that the distribution of r|^^| actually 
has an exponential tail. 

The functions are central to the argument to come. First, we show 
that they are bounded and Lipschitz, with appropriate constants. 



h-diff 



Lemma 2.4 For all j € C, 

\h'}{j) - h>}{s)\ < 2\\f\\T/p. 
Proof. For any j G C, we can write 

/•oo 

-h'^ij) = E,{(/(Xn-vr'^(/))/[rf,^<t]}cit 

/oo 
E,{(/(Xf)-7r^(/))/[rf;j>t]}dt. 

Then, by the strong Markov property, we have 

E,{(/(Xn-^^(/))/[r{;j<t]} = E,|E,{(/(Xf_^,^)-vr'^(/))/[Tf^^<t]}| 



g%{v)¥.,{f{Xl^)-^^{f))dv, 



where g'^g denotes the probability density of the random variable t'^^^ for the 
process started at j. Hence it follows that 

roo 

E,{(/(Xf)-^^(/))/[r|;^<t]}dt 



oo fOO 



dt / dvg%{v)^,{f{Xl^)-7:^{f))l{v < t}. 
JO 



Now, since 



l{v<t}\E,{f{Xt_,)-7r'^{f))\dt 

/•oo 

< 211/11 / dM^{XnX^^{0) = s),7Tn < oo, 
Jo 



we can use Fubini's theorem to conclude that 



E,{ifiXn-7r^{f))I[Tf,y<t]}dt 



g^^{v)\ I E,{fiXU-^''if))dt\ dv 



g%{v^{s)dv = h'}{s). 



Hence, from ()2.6p and Lemma 12.21 it fohows that 

\h1ij) - h'^is)] < 211/11 E,-(t{;P < 2\\f\\T/p, 

as required. 



h-equation 



Stein-eq 



In particular, the function is itself bounded. 



A similar argument, by conditioning on the time of the first jump, shows 
that 

(2.7) h'^U) = -q-'{f{j)-n'^{f)}+ J] qj\^kh%^), 

where qj := Ylk£Cu{o} ^jk < oo because X is conservative, and the sum 
in (|2.7p is absolutely convergent because /i^ is bounded. This can be rewrit- 
ten in the form 

(2.8) {Q'^h^fm = f{j)-^^{f), j^C, 

so that, for any bounded / and for any probability measures n and u on C, 
we have 

(2.9) vr^(Q"/i^) = TT'^in-TT^f). 

In the terminology of Stein's method, (j2.8p determines /ij to be the solu- 
tion h of the Stein equation {Q'^h){j) = /(j)— 7r^(/) for the distribution tt^, 
corresponding to the given function /. Also, by Dynkin's formula, we have 



Dynkin (2.10) 



^'^(Q^/i) 
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for any bounded function h (for the special case h = hp this follows 
from ()2.9p ). These considerations put us into a position to prove Theo- 
rem [27TJ 

Proof of Theorem 12.11 Take any probability measures /i and ly on C. 
Then (gS]) gives 

whereas (|2.10p gives TT^^{Q^h'j) = 0. Taking the difference, we obtain 
Pi-diff| (2.11) vr'^(/)_vr-(/) = n^^iQ'^h'} - Q^^^). 

Now, for bounded h and for any i €z C, 

{Q^h-Q^^h){i) = Y,qUKk)-h{i))-Y,q^Mk)-hii)). 

keC k&C 

with both sums absolutely convergent, and, from (j2.ip . it then follows that 

(Q^/i-Q^/i)(i) = qioY.{u{k) - ii{k)){h{k) - h{i)) 

kec 

= qioY.{u{k)-fi{k)){h{k)-h{s)), 
kec 

since X^feeC ^^i^) ~ X^fcec l^i^) ~ ^- Hence, from (|2.11j) . we have 

T^'iD-^-'if) = ^vr'^(.)g.oEM^)-MA:))(/iK^)-^/(^))' 

iec kec 

and, from Lemma 12.41 this gives 

penultimate I (2.12) |vr^(/) _ vr''(/)| < ^TT>'{i)qio2\\f\\{T/p)\\iy - ij\\tv 

iec 

Thus it follows that 

(2.13) Wtt" - nf^WTv < i2T/p)^7r^{i)qio\\iy- fiWrv, 

iec 

and (|2.4p then implies that 

(2.14) IItt'" - Tr'^llTy < i2UT/p)\\iy - hWtv- 

This, by the Banach fixed point theorem, establishes the first part of the 
theorem, and the second follows by taking u = m, and using the fact that. 



contraction 



final 



for probability measures F and G, cZtv(-P) G) = ^'^F — G\ 



TV- 



quasi-longtime 



We now turn our attention to the distribution of X{t) for fixed values 
of t, starting from any initial distribution, and compare it to m, the distri- 
bution at any time of the return process X™ started in the quasi-stationary 
distribution m. We begin by taking the initial state of X to be s, and re- 
mark later that this restriction makes little difference, provided that s is hit 
at least once. 

Theorem 2.5 Let B := Tqg/p > 1. Then, under Condition A and if also 
2UT/p < 1, there is a constant K such that 



d^y{£s{Xit)),m) < Ut + KBj- + {2/e)P^/^^^ =: r]{t). 

pt 



Remark. Hence, if UB^T/p <C 1, the distribution is close to m 

for all times t such that 

B'^T/p < i < f/-\ 

Proof. The argument is based on coupling two copies X^^^ and X(^) of 
the return process X™, with X^^^ in equilibrium and with X^"^^ starting 
in ,s. The coupling is achieved by forcing X^^^ to follow the same sequence 
of states as X^"^^ after the first time that it hits s, and to have identical 
residence times in all states other than s. Define 



rlyil) := inf{t > 0: X^'\t) = s}; r°,}(2) := 0, 



and let 



rf,}(/) := inf{t > r^-^l): ^ X^'\t) = s}, 1 = 1,2, 

denote the n-th return time of X^''^ to s. Then, conditional on the event 
that T"{'5}(1) = we have 

d,v(/:(Tf,}(l) I T^^il) = «),£(Tf,j(2))) 

= dTv{Sy*q~^G{n,l),qj'^G{n,l)) < caqsvn'^^'^, 

for a suitable constant cq, where G{n, 1) is the Gamma distribution with 
shape parameter n and unit scale parameter, and dy is the point mass at v. 
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Ist-bnd 



gamma- 


-bnd 








2nd- 


-bnd 








3rd- 


-bnd 



Hence, for any n > 1, we can couple X^^^ and X^'^^ by arranging that 
T^^y{l) = T^^y{2), with the two processes to be run identically thereafter, 
and the probability of this coupling failing, conditional on T|^|(1) = v, is at 
most CGqsVn~^^'^ . Thus, in particular, 

(2.15) d^^{C{X^'\t)),C{X'^'\t))) < F[Tf^y{2)>t]+CGqsn-'/^T/p, 



using Lemma 12.21 It now remains to show that we can reach the bound 
given in the theorem by choosing n almost as a multiple of t. 

Now t"^|(2) is a sum of independent random variables, each with distri- 
bution £s(T|^|), where t"^, is defined as in Lemma |2.2[ By that lemma and 
Markov's inequality, it follows that 

Fr[T^^^>2T/p] < 1/2, reC, 

and hence that 
¥,[T?l.>2kT/p] < 2-^ 



exp 



plog2 f2kT\ 
2T V W 



for all /c > 1. 



Thus the distribution Cs{T^y) is stochastically bounded above by that of 



2T 

P 



1 + 



log 2 



-E 



where E has a standard exponential distribution. Hence the distribution 
of rj^^|(2) is stochastically bounded above by that of 



2T 

V 



n + 



1 



log 2 



Gr 



where Gn ~ G{n, 1). The inequality P[G„ > 2n] < (2/e)" thus implies that 



r{"}(2) > 



2Tn 



1 + 



log 2 



< (2/e)^ 



So, for any fixed t, using 1 + 2/ log 2 < 4, we take n 
in (j2.15p . giving 



nt := [tp/8T\ 



(2.17) d^^{C{X^'Ht)),C{X^'Ht))) < {2/er+CGqsn;'^'T/p, 



from which it follows that, for t > 16T/p, 

(2.18) d^y{C{X(^\t)),C{X^^\t))) < {2/er/'^^ +4cG 



{Tqs/pf'\ 

Vq^ 
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We first observe that C{X^-^\t)) = m for all t. Then we have 

P[r{o}(l) <t] = < Ut, 

where r{o}(l) := mf{t > 0: X'^'^^t) = 0}. On the event that X^^^ and X^^) 
are successfully coupled at t"J| < t, it thus follows that the event that nei- 
ther hits before t has probability at least 1 — Ut, and, on this event, X^'^\t) 
is also the value of an X-process starting in s, since has had no visits 
to before t. This, together with ()2.18p . completes the proof. I 



Remark. Denoting by j4({s}, {0}) the event that X hits s before 0, the same 
argument can be used to show that dTy{£-kiX{t) \ A{{s},{0})), Cs{X"^{t))) 
is at most ri{t) for any k G C, under the conditions of Theorem 12. 51 Hence, 
conditional on the event that X hits s before reaching 0, the distribution 
of X{t) starting from any k £ C is also close to m for all times t such that 

B^T/p ^ t ^ U-\ 

provided that UB'^T/p < 1. Thus the quasi-stationary distribution m is 
then indeed the appropriate long time approximation to the distribution 
of X in C, for times t < U^^. 

Note also that the coupling used in Theorem 1 2 . 5 1 may be very pessimistic, 
only making use of the residence times in s. For most processes, the variabil- 
ity in the remaining residence times and in the possible sequences of states 
can be exploited to get sharper bounds. However, in the examples for which 
we make computations below, the quantity B'^T/p is of only polynomial 
order in the size of the system, whereas is exponentially large; hence 
even this crude estimate is more than adequate. 

3 Birth and death processes 

bfed| 

Consider now a birth and death process with C = {1, 2, . . . , N} (for C = N, 
replace A'" by oo in what follows) having birth rates bj > 0, 1 < j < N, 
with 6o = and 6jv = 0ifA^<oo, and with strictly positive death rates 
dj, j G C. It is convenient to introduce the quantities {aj, j G C), where 
ai = 1 and, for j > 1, 

6i • • • bj_i 
d2---dj 
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The return process with /i = S^i}, equivalent to re-defining di to be zero, 
is then recurrent if a+ := Ylj>i'^j < ™ which case vr^(j) = aj/a+, 
so that its computation is very easy. We now wish to investigate when 
this distribution can be used as a reasonable approximation to the effective 
steady state behaviour of the process. 

In order to apply Theorems 12. II and 12. 51 we need to choose a state s > 1, 
and find values for p,T,B and U. For p, let r^, A; > 1, be the probability 
that the process starting in k hits s before it hits 0, where s > 1. If A; > s, 
then rfc = 1. Otherwise, rg = 0, = 1 and 



{bk + dk)rk = hkTk+i + dkVk-i, 
leading to = ak/cTs, where 

k ^ 

0-0 = and ak = ^ — 

Since ak is non-decreasing in k, we can take 
p-bdl (3.1) p = n = l/{dias), 



l,2,...,s-l. 



for A; = 1, . . . , s. 



for any state s £ C. 
For T, we note that 



s-l 



Tl-def 



(3.2) Ti := rnax Efc{r|,^o}} < YjTT^^j < 

and that, for k > s, 

k fc ^ oo 



j=s+l 



dnUj 



ai 



(Anderson 1991, Chapter 8). Since the latter quantity is increasing in k, we 
may take 



pkp2 1 (3.3) T := max(Ti,T2), where T2 := y^aj. 



d^aj 



Note that then Condition A(ii) holds if T2, the so-called "D series", con- 
verges, and that T2 < 00 is a necessary and sufficient condition for a 
birth and death process to have a unique quasi-stationary distribution (van 
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Doorn 1991, part 2 of Theorem 3.2). Note also that B := Tqs/p can be 
bounded using (j3.ip - (|3.3p , together with the fact that Qs = bg + dg. 
Finally, the quantity U can be evaluated as 

(3.4) [/ = 



(6i + di) Ei(r{i,o}) 1 + h E2(r{i}) ^f^^ a, ' 
because, also from Anderson (1991, Chapter 8), 

1 ^ 1 ^ 

and in particular, since ai = 1, 

N N 

1 + 6iE2(t|i|) = 1 + J]a, = J] a,. 

i=2 i=i 

In order to apply Theorems 12.11 and 12.51 in practice, we need to be able 
to bound the quantities p, T, B and U by assigning concrete expressions in 
terms of the hj and dj to replace (j3.ip - (l3.4p . Simple estimates can be derived 
under the assumptions that the death rates dj are increasing in j, and that 
the ratios hj/dj are decreasing, with hi/ di > \. If this is the case, define 
s > 1 in such a way that hs/dg > 1 > fts+i/c^s+i) and let 1 < si < s < S2 be 
such that 

hs^/ds^ =: pi > 1 > P2 := hs^/ds^- 
Then Xj := dj^iaj+i/di is maximal at j = s, and 

Xj > pI\ Sl<j<S; Xj/Xs < /02~''^ i > •S2- 

Hence, from (I3.1D-(I3.4D. we have the bounds 



s-l 

p 



u < dsM-^)pI''; 

T2 < {s2-s + i/{i-p2)} ^; 

j>s+l ^ 

s—1 ^ i s — 1 ^ 

j=l ' ' j=0 i=l * 
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Thus if, for instance, dj grows at most polynomially fast in j, with the 
sum ^j>i dj^ < oo, and if si and s — si are large and of comparable size, 
then T/p is roughly of order and B = Tqg/p of polynomial order in s, 
whereas U is geometrically small with s, making UT/p very small indeed. 

More precise calculations for the stochastic logistic model of (|l.ip as 
A — )• cxD, with s = [kA\ , give 

T = 0(log A), p > l-d/b, U < h+ ^1 , B = O(^log^), 

so that UT/p is geometrically small in A as ^ — t- oo. Thus, for the stochastic 
logistic model, the unique quasi-stationary distribution can be very closely 
approximated by any return distribution, as long as A is large. Entirely 
similar estimates are true for the SIS epidemic model, which models the 
number of susceptibles in a closed population of size A^, to be thought of as 
large but finite. The process is a birth and death process on {0,1,..., N} 
having rates 

(3.5) hi := \i{l-i/N) and di := i^i, 0<i<N; 

in this case, UT/p is geometrically small in if < A, and there is a 
quasi-stationary distribution close to s := [A^(l — /i/A)J. 
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